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Abstract. There has been considerable interest recently in the generation of 
azimuthal phase functions associated with photon orbital angular momentum (OAM) 
for high-dimensional quantum key distribution (QKD). The generation of secure 
quantum keys requires not only this pure phase basis, but also additional bases 
comprised of orthonormal superposition states formed from the pure states. These 
bases are also known as mutually unbiased bases (MUBs) and include quantum states 
whose wave functions are modulated in both phase and amplitude. While modulo 
2ir optical path control with high-resolution spatial light modulators (SLMs) is well 
suited to creating the azimuthal phases associated with the pure states, it does not 
introduce the amplitude modulation associated with the MUB superposition states. 
Using computer-generated holography (CGH) with the Leith-Upatnieks approach to 
hologram recording however, both phase and amplitude modulation can be achieved. 
This paper presents a description of the OAM states of a 3-dimensional MUB system 
and analyzes the construction of these states via CGH with a phase modulating 
SLM. The effects of phase holography artifacts on quantum-state generation are 
quantified and a prescription for avoiding these artifacts by preconditioning the 
hologram function is presented. Practical effects associated with spatially isolating 
the first-order diffracted field are also quantified and a demonstration utilizing a liquid 
crystal SLM is presented. 
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1. Introduction 

With advancements in digital electronics, the principles of holography [1] continue to 
be utilized for new approaches to imaging and diffractive control of optical fields. In 
digital holography (DH), an imaging sensor and computer record the hologram and 
image reconstruction is accomplished via digital processing techniques to present high 
resolution digital images. [2] In computer-generated holography (CGH), a hologram 
function can be generated computationally and then displayed via digital printing, 
photolithography or display-based technologies. [1, 3] A 2006 special issue of Applied 
Optics presents an overview of recent advances in DH and CGH as applied to microscopy, 
optical trapping and tweezing, interferometry, data recording, beam shaping, and 
metrology. [4] 

There has been considerable interest recently in the generation of azimuthal 
phase functions associated with photon orbital angular momentum (OAM) for high- 
dimensional quantum key distribution (QKD).[5] The complex fields associated with 
these functions are of the form exp(U9) where 9 is the azimuthal coordinate. The 
integer / denotes the photon orbital angular momentum in units of Plancks constant 
over the range — oo < / < oo. While the conventional BB84 protocol for QKD [6] is 
based upon photon spin angular momentum and utilizes two polarization bases (e.g. 
the horizontal and diagonal bases), extensions of this protocol to higher dimensions 
is possible. [7] Whether in two dimensions or higher dimensions, secure transmission of 
quantum keys is ensured by the use of mutually unbiased bases (MUBs).[8, 9] The states 
associated with each MUB are equally weighted superpositions of the states associated 
with the other bases. The measurement of a single photons state with an incorrect choice 
of basis results in an unbiased, or meaningless, probabilistic outcome. Preserving this 
essential feature of the BB84 protocol while implementing QKD with higher dimensional 
OAM Hilbert spaces, requires one to generate and discriminate each state in each of the 
MUBs associated with OAM. These OAM MUB states correspond to optical fields that 
are modulated in both phase and amplitude. 

Recent demonstrations5 of the generation of photon OAM have utilized high- 
resolution programmable optical path modulators to impart the azimuthal phase 
dependence of the pure states onto optical wavefronts. For values of / > 1, the limited 
range of these spatial light modulators (SLMs) requires that the optical phase function 
be introduced modulo 2n. Assuming monochromatic conditions and adhering to the 
sign conventions of Goodman [10] wherein a positive optical phase delay introduces a 
negative phase error, the complex transmittance associated with an optical phase delay 
function 16 introduced modulo 2n can be shown to be [11]- [13] 

oo 

t(0)= sinc(l-m)e imW (1) 

m=— oo 

where m is an integer and the sine function is used according to the convention 
sinc(x) = sin(irx)/iTx. The amplitude diffraction efficiency, given by the sine function, 
is spatially uniform and unity for m = 1 and zero otherwise. While modulo 2n optical 



Holographic generation of complex fields 



3 



phase control is well suited to creating the azimuthal phases associated with the pure 
states with high diffraction efficiency, the diffractive optics function is calculated without 
the benefit of a reference wave to encode amplitude information. 

Methods for computing hologram functions that give rise to complex field 
modulation include both analytic and numerical techniques. Historically, the simplicity 
of fabricating binary holograms motivated the development of numerical techniques for 
calculating and optimizing the binary hologram function. [1, 14] Recent developments 
in SLM technologies allow analytic hologram functions to be displayed in real time 
with megapixel spatial resolution and 8-bit phase resolution. Consequently, analytic 
expressions for hologram recording functions can be used directly in CGH. The off- 
axis hologram recording approach introduced by Emmett Leith and Juris Upatnieks 
[15, 16, 17] defines a useful approach to CGH with SLM technologies providing both 
phase and amplitude modulation of the diffracted field and allowing separation of 
the diffracted orders. This approach also allows the hologram transmittance and the 
complex fields associated with the diffracted orders to be expressed analytically for 
purpose of analysis. 

While photon polarization supports QKD in a 2-dimensional Hilbert space, in 
this paper we present a description of OAM MUBs for a three-dimensional Hilbert 
space and consider the implications of generating these quantum states via CGH with 
SLM technology. While our analysis is valid for higher dimensional QKD, we present 
this next-higher dimensional generalization for clarity. Toward this end, the complex 
transmittance associated with the Leith-Upatnieks approach to recording thin phase 
holograms is reviewed and evaluated for the case of holographic generation of a MUB 
state. Artifacts in the holographically generated complex field are discussed and their 
effects on a QKD system quantified by calculating the weighted inner product of the 
holographically generated field with the theoretical MUB field. Through preconditioning 
the CGH function, these artifacts are eliminated. Finally, a demonstration is presented 
in which an extended graphics array (XGA) format liquid-crystal-on-silicon (LCOS) 
SLM is utilized to holographically generate an optical field associated with a three- 
dimensional OAM MUB state. 

2. Mutually Unbiased Bases of OAM States in 3 Dimensions 

A 3-dimensional Hilbert space admits a maximum of four MUBs. Each of these bases 
is, by definition, orthonormal and contains three basis vectors. We may freely choose 
as one of these bases, designated MUBq, the three pure OAM states corresponding to 
an angular momentum, l a = ah, lb = bh, and l c = ch, where a,b,z G Z. To further 
reinforce the vector nature of these states, we utilize Diracs ket notation [18] to provide 
an abstract description of each of these three basis vectors in this MUB. 

MUB = {\a),\b),\c)} (2) 

These three ket vectors are orthogonal, providing a ^ b ^ c , and span the 3-dimensional 
Hilbert space. Any other quantum state \ip) in this space can be written as a linear 
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b) 
c) 




ai) oc 
61) oc 
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a) + |6) + |c) 
a) + z&) + z 2 
a) + ^ 2 |6) + z 


c) 
c) 




a 2 ) oc 
b 2 ) oc 
c 2 ) oc 


a) + |6) + z 
a) + z\b) + 
a) + z 2 \b) + 


c) 
c) 

z 2 \c) 




03) oc 
63) oc 
c 3 ) oc 


a) + |6) + z 2 
a) + z\b) + z 
a) + z 2 \b) + 


c) 
c) 
c) 



Table 1. Dirac notation description of mutually unbiased bases in 3-dimensional 
Hilbert space. 



superposition of the three basis vectors, 

|^) = a \a) +P\b) +t|c), (3) 
with a, 7 G C and 

a*a + + 7*7 = 1 (4) 
We require that these states be normalized with inner product 

{Aj) = kj for h -3 e {a, 6, c} (5) 
Consider a second MUB basis in this 3-dimensional Hilbert space, MUB 1 = 
{|ai), |ci)}. Two characteristic features are needed to define another MUB basis. 
First, the three vectors in such a basis must be orthonormal and second, any state in 
MUBi should be equally distributed over the three vectors in MUB . Hence, 

l(%')| 2 = ^ fori, j G {ai, b±, Ci} (6) 

The only four MUB bases in this 3-dimensional Hilbert space are shown in Table 1. 
Note that in Table 1, z — exp(i2/3) and we have suppressed the normalizing factor of 
1/V3 in each of the nine MUB states in the last three columns. 

Following the Dirac notation, we can provide a configuration-space representation of 
each of these twelve azimuthal MUB states in terms of a real- valued amplitude function 
A{6) and a real-valued phase function The MUB states are purely complex and 

may be expressed as phase function in the optical plane. 

(0\a) = e iae (7) 

(9\b) = e M (8) 

(9\c) = e icd (9) 

Each of the other nine states in Table 1 is of the form, 

|d)=a|a)+/3|6)+7|c) (10) 

and accordingly, 

(9\d) = A(9)e~ im (11) 

The corresponding real- valued amplitude and phase functions can be expressed in terms 
of Eq. 11 as 



A(9) = ^(d\d) (12) 

m = -tan (13) 
1 ; \Re((9\d))J 1 ; 
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For the purpose of this paper, and without loss of generalization, it is convenient 
to assign the pure states \a), \b) and |c), the quantum numbers |1), |0) and | — 1), 
respectively. Furthermore, the analysis that follows will concentrate on CGH generation 
of one representative MUB state, namely 

| C 3) = ^(|l) + e-^ 3 |0) + |-l)) (14) 

which may be expressed in complex field notation with the following real-valued 
amplitude and phase, 



A{6) = -=Vl -2cosfl + 4cos 2 fl (15) 
v3 



$(0) = - tan" 1 



V3 



4 cos (6)-lJ 

3. Complex Field Modulation via Phase Modulation CGH 



(16) 



The off-axis hologram recording approach introduced by Leith and Upatnieks defines a 
useful approach to complex field modulation by CGH with SLM technology. Introducing 
the reference wavefront off axis leads to spatial separation of the diffracted orders 
allowing the first-order diffracted field to be isolated from other diffracted components. 
The functional dependence of the complex field associated with each diffracted order 
may be derived explicitly. Fourier series representations of thin phase hologram 
transmittance are given in the standard references for the case of one-dimensional 
modulation introduced by two planar wavefronts.[l, 10] In this section, we present the 
more general form that includes spatial modulation of the recorded field and show that 
phase modulation holography leads to artifacts in both the amplitude and phase of the 
holographically generated optical field. Based on analytic expressions for these artifacts, 
we present a prescription for preconditioning the computer-generated hologram in order 
to avoid these artifacts in the diffracted field. The wavelength dependences associated 
with holography are suppressed by assuming monochromatic conditions throughout. 

3.1. Two- Dimensional Hologram Recording Function 

Let a(p, 6) and <p(p, 6) represent spatially varying real-valued amplitude and phase 
functions, respectively, where p and 9 are the transverse radial and azimuthal 
coordinates, respectively. The complex optical field comprising both a(p, 8) and 4>(p, 6) 
may be written as 

E(p,0) =a(p,#)e W) . (17) 

Similarly, the phase associated with an off-axis reference wavefront of slope a along the 
Cartesian y direction is k ay which, in polar coordinates can be written as kp sin 9, 
where k is the wave number. Hence, the reference wave R(p, 6) is written as, 

R(p,9) = be- ikapsine , (18) 
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where b is a constant amplitude. In order to maximize fringe contrast and maintain 
a monotonic relationship between fringe contrast and amplitude, the reference beam 
amplitude b is chosen to be equal to the maximum value of the amplitude a(p, 9) denoted 
Umax- The spatially varying intensity associated with the superposition of both fields is 

I(P, 8) = a 2 (p, 9) + a 2 max + 2a(p, 9)a max cos (0(p, 9) - hap sin 9) (19) 

. To allow for scaling of the hologram recording function, the computer-generated 
holography function f cg h(p, 9) is defined as 

/c 9 .M)=<4^, (20) 

1-max 

where I max normalizes I(p, 9) to unity and the parameter a establishes the maximum 
value. 

In principle, this function can be represented as a thin amplitude hologram via an 
amplitude modulating SLM or as a thin phase hologram via an optical path modulating 
SLM. Commercially available optical path modulating SLMs include piston-only micro- 
mirror arrays [19] based on micro electro-mechanical systems (MEMS) manufacturing 
and parallel-aligned nematic liquid-crystal (PA-NLC) media when used with linearly 
polarized light aligned with the NLC director. [20]- [24] The variable birefringence of PA- 
NLC media also allows them to be used in conjunction with orthogonal polarizers to 
create amplitude modulation and is utilized in this manner for display applications. It 
should be noted however that the complex transmittance also includes a parasitic phase 
modulation term. Marquez et al. [25] have shown that twisted NLC (TNLC) media 
can be configured with a tandem arrangement of wave plates and polarizers to produce 
either nearly pure amplitude modulation or nearly pure phase modulation. 

3.2. CGH with an Optical Path Modulating SLM 

For the case of implementation with an optical path modulating SLM, the specific form 
of the complex transmittance is determined by the phenomenology of the SLM. In Sec. 5, 
we will sample and scale the CGH function of Eq. 20 in order to define an 8-bit gray 
scale XGA signal for the SLM driver. The local response of the LC media is such that 
the refractive index, and correspondingly the optical phase retardance, decreases with 
increasing intensity as defined by Eq. 20. Therefore, the change in phase retardance 
introduced by the SLM is, 

A^ S lm(p,9) = -a 1 -^- (21) 

-Lmax 



where a defines the maximum change in phase introduced by the SLM. The resulting 
complex transmittance t(p, 9) of the hologram is given by, 

t( p ,9) = e' la ^. (22) 

By substituting Eq. 19 into Eq. 22 and applying the Jacobi- Anger formula along with 
the identity30 

Jm{-Z) = {-l) m J m {z), 
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the transmittance can be written in terms of the diffracted orders as, 

m=~oo \ a max ) 

where J rn is the Bessel function of the first kind of order m with the maximum value 
of the argument given by a' = 2aa 2 nax / I max . The coefficient to the series is a spatially 
uniform phase shift. Assuming that the SLM is illuminated with the reference wave of 
Eq. 18 with amplitude b now taken to be unity, the transmitted optical field is given by 
the product of Eqs. 18 and 23, 

E out (p, 9) = £ (-1)- Jm L> a -^\ e -«M^), (24) 

m=-oo \ a max ) 

where the optical phase of the m-th order diffracted component is 

m (p, 9) = -m(j>(p : 9) + (m+ l)kap sin 9 + °- (25) 

Z y d max J 

It is evident that the recorded phase <p(p, 9) is found in the m — — 1 order. Applying 
the Bessel function identity J m (z) = (— l) m J m (z),[26] the spatially varying amplitude 
of the m — — 1 diffracted order may be written as, 

a_ 1 (p,9) = J 1 (a' a -^l) (26) 

\ ^max / 

The diffracted amplitude ai(p,9) is related to the recorded amplitude a(p,9) through 
the J\ Bessel function leading to nonlinearities in the amplitude of the holographically 
generated optical field. The spatially varying efficiency with which optical power is 
diffracted into this order is given by the square of a\(p, 9), 

,. l(p ,e) = Jl , (27) 

and achieves a maximum value of 33.9% at the peak of the J\ Bessel function. If a 
given value for the amplitude of the diffracted field is to be defined by a unique value of 
the argument of the Bessel function, then the parameter a' should be chosen such that 
the Bessel function remains monotonic, a' < 1.84. The upper limit for the maximum 
intensity I max is ^a? max . For this case, a' = a/ 2 and the Bessel function will remain 
monotonic if a < 3.68. 

The spatially varying phase associated with the m — — 1 diffracted order, (f>i(p,9), 
is given by 

<P- l (x,y)=<P(p,9) + ^(^±) 2 (28) 

While Eq. 28 indicates the desired phase <f>(p, 9) is constructed by the hologram, it 
also shows the presence of the aberration term %r ( a ^ p ' 9 ^ ) . This aberration term is 
proportional to the spatially varying intensity a 2 (p,9) and assumes its maximum value 
when a(p,9) = amax. For a' = 1.84, the maximum phase error would be 0.92 radians 
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or 0.15 waves, within the classical quarter-wave Rayleigh criteria for diffraction limited 
performance. [27] 

For applications requiring a high degree of fidelity in the holographically generated 
wavefronts, CGH offers the opportunity to precondition the computed hologram to 
compensate for both the J\ nonlinearity and the phase aberration. In order to generate 
an optical field with modulated amplitude A(p, 9) it is only necessary to define the 
CGH amplitude function a(p, 9) such that the amplitude of the m — — 1 diffracted field 
satisfies the relationship 

o_i(p, 6) = J\ L ,a ^\ = cA(p, 9) (29) 

\ Q>max J 

where c is a scaling factor that for a given value of a' matches the maximum values of 
Ji and A(p,0). Given the preconditioned amplitude function a(p,9), it is then possible 
to define the preconditioned phase function <j)(p,9) that gives rise to 9). This is 
done by setting <p\{p, 9) in Eq. 28 equal to $(p, 9). 

<Mx, y ) = 0M) + ^(^^) =*(p,e) (30) 

z \ a max j 

4. Generation of an OAM MUB State via Phase Modulation Holography 

The MUB states that we generate using CGH will never exactly represent the theoretical 
MUB state. In a QKD system, numerous effects will introduce cumulative deviations 
from the theoretical state. In addition to the artifacts of phase modulation holography 
described above, these can include optical aberrations within the system including 
atmospheric aberrations. [28] In this section, we evaluate the effects of phase modulation 
holography artifacts on the fidelity of a three dimensional MUB state. We also consider 
the case where the hologram is preconditioned to minimize these artifacts and evaluate 
additional diffraction effects associated with spatially isolating the first diffracted order. 



4-1. Artifacts of phase modulation holography 

For the purpose of illustration, consider the \cs) state of the previous section. Recall 
from Eqs. 15 of Sec. 2 that the amplitude and phase associated with the superposition 
state may be written as, 

A{9) = -i= (l -2cos(#) +4cos 2 (fl))^ (31) 

$(#) = - tan ( ^ ) (32) 

These are shown in gray scale in the first row of Fig. 1 with the amplitude plot normalized 
to unity and the phase plot in units of radians. Taking the amplitude a(p, 9) and phase 
4>(p,9) in the CGH calculation to be exactly those of the |c 3 ) state, a(p,9) = A(p,9), 
amax = Amax, and <j)(p, 9) = $(p, 9), we can evaluate the effect of the phase modulation 
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holography artifacts on the holographically generated |c 3 ) state. Eqns. [22, 24] give the 
amplitude and phase of the m — — 1 order diffracted field to be, 

... P3) 

<t>-m=m^{^-) 2 (34) 

Upon substitution from Eqs. 33 the amplitude and phase of the holographically 
generated field are, 

a_i = Ji ( (l - 2 cos {&) + 4 cos 2 (0)) 3 ) (35) 

*->=- tm {idhi) + ^ (1 " 2 cos (9) + 4 cos2 (36) 

The amplitude and phase functions of Eqs. 35 are evaluated for maximum values of 
the Bessel function argument given by a' = 0.610, 1.84, 3.13, and 3.83 corresponding 
to the first half-maximum, peak, second half-maximum, and second zero of the Bessel 
function, respectively. The results are shown in gray scale in the remaining rows of 
Fig. 1. For the case of a' = 0.610, the Bessel function is nearly linear and the first-order 
amplitude is nearly identical to the \cs) state amplitude shown in the first row. For 
a' = 1.84, the nonlinear scaling is more significant and the effects become discernable 
in the amplitude plot. For larger values of a', the non-monotonic nature of a\(Q) is 
clearly evident. At these larger values of a', the phase aberration is also discernable in 
the plots. Quantum mechanics provides us with a physically-meaningful metric over the 
Hilbert space spanned by the apertured OAM quantum states |Z), where I e Z , that 
allows us to compare the fidelity, or the angle between ket vectors, of our holographically 
generated optical MUB states with the theoretical states. In particular, the metric takes 
the form of an inner product over the ket vectors between our holographically produced 
MUB state \ip h ) and the theoretical MUB state \ipt), each of which can be represented 
as wave functions in the azimuthal plane {p, 9}. 

mh) = mh\ (37) 

Here, a is the aperture radius, and is the norm of the wave function similarly defined 

as, 




2tt 

ifj*{p,6)ij{p,6)pdpd6 (38) 

10 jo 

This inner product in the Hilbert space represents the cosine of the angle in Hilbert 
space between our observed and idealized MUB state but, more importantly, the square 
of this inner product represents the probability P that our holographically produced 
MUB state will be measured (by a perfect MUB state detector) in the theoretical state. 

P=\&t\ii>h)\ 2 (39) 
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pupil plane pupil plane 

amplitude phase 




Figure 1. Calculated amplitudes and phases comparing theoretical complex field to 
that generated by CGH with phase modulation holography for various values of the 
phase scaling parameter a'. 



We will examine the dependence of this fidelity measure on 1) the phase parameter a' 
for the theoretically defined first diffracted order and 2) the reference wavefront slope 
for numerical simulations of generating and spatially isolating the first diffracted order. 

Fig. 2 shows the calculated probability function as a function of the parameter a'. 
For sufficiently large values of the tilt a, the maximum intensity I max can be made 
arbitrarily close to 4A^ a;E . For the purpose of Fig. 2, it is sufficient to choose a to 
correspond to 100 waves of tilt across the aperture and let the maximum intensity I max 
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1 2 3 

phase scaling, a' 



Figure 2. Calculated inner-product-based probability quantifying the fidelity with 
which the complex field is generated via CGH vs. the phase scaling parameter a' for 
each of the mutually unbiased bases in 3-dimensional Hilbert space. 



be represented by 4a^ a;I ,. The maximum value for the argument of the Bessel function 
then becomes a' = a/2 and the peak of J\ occurs when A(9) = A max and a' = 1.84 
or, equivalently, a = 3.68. The calculated probabilities for all states within a given 
MUB are identical. Since the pure states of MUB have constant amplitude, the phase 
aberration term of Eqs. 33 results in a constant phase shift. Similarly, the 3\ Bessel 
function is evaluated at a single value of the argument thus avoiding nonlinearities 
in the amplitude response. Consequently, the calculated probabilities for MUB are 
unity for all values of a'. The probabilities associated with the superposition states of 
the remaining bases all decline with increasing a' as shown with the states of MUB2 
and MUB3 falling on the same curve. At a' = 1.84, the maximum value for which 
Ji remains monotonic, the probability functions have decreased to about 0.90. As a' 
increases through the range where Ji is no longer monotonic to a' = 3.83, the second 
zero of Ji, the probability functions decrease further to the vicinity of 0.3 to 0.5 with 
the states in M\JB X experiencing the largest reduction. 



4-2. Effects of Spatially Isolating the First- Order Diffracted Field 

Another practical effect associated with the holographic generation of optical fields is 
that of isolating the first-order diffracted field from other diffracted orders. In practice, 
this is accomplished by choosing the angle of the reference wavefront to be sufficiently 
large that the angular diffraction is greater than the divergence associated with the 
diffracted fields. The first-order component may then be transmitted by an aperture 
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while blocking the other components. 

The effects of spatially isolating the holographically generated MUB 3 field may 
also be evaluated numerically. The first row of Fig. 3 again shows the amplitude and 
phase associated with the theoretical function as defined in Eqs. 31. Note the gray scale 
for the phase plots is different than that in Fig. 1 in order to remain consistent with the 
results that follow. The second row shows the numerically calculated holographically 
generated optical field for the case of a reference wave with 10- waves of tilt. In this case, 
the amplitude and phase artifacts of phase holography have been avoided by using the 
preconditioned amplitude and phase functions a(p, 6) and 0(p, 6) as defined in Eqs. 29,30 
to calculate the CGH function of Eq. 20. The parameter a' is taken to be 1.72 in order 
that the Bessel function may be approximated by a polynomial and Eqs. 29,30 may be 
solved for a(p, 6) and <p(p, 8), respectively. The complex transmittance of Eq. 22 is then 
used to define the optical field transmitted by the hologram as elements of a 768 x 1024 
numerical array. The far-field amplitude is calculated by FFT and is shown in the 
first column of Fig. 3. The gray scaling has been chosen to be nonlinear to emphasize 
the m — ±1 diffracted orders in the figure. A rectangular aperture centered at the 
m — — 1 diffracted order is chosen such that the width of the aperture is equal to the 
order spacing. The field transmitted by the aperture is then transformed by FFT back 
to the pupil plane. The resulting amplitude and phase are also shown in the second 
and third column, respectively. Note that while the orders appear to be well resolved, 
artifacts occur in the amplitude and phase functions. These artifacts may be attributed 
to interference effects from adjacent orders and diffraction due to the finite aperture 
size. The third and fourth row show the effects of increasing the reference wave tilt to 
50 and 100 waves respectively and scaling the aperture size accordingly. Note that the 
artifacts are significantly reduced. 

The fidelity of the holographically generated fields may again be quantified by 
calculating the probability function in Eq. 39 from numerical data similar to that shown 
in Fig. 3. The results are shown in Fig. 4 for the case of the \c^) state and are calculated 
for the cases of 10, 30, 50, 75, and 100 waves of reference wave tilt. The squares 
represent the case where the CGH function was not preconditioned and therefore include 
the effects of the amplitude and phase artifacts described in Figs. 1,2. The diamonds 
represent the case in which the artifacts are precompensated. Fig. 4 shows that the 
combined effect of preconditioning the hologram and increasing the reference tilt from 
10 to 100 waves is to increase the inner-product-based probability from 88% to 99.6%. 

5. Experimental Demonstration with a LCOS SLM 

In this section, we utilize a LCOS SLM to demonstrate holographic generation of 
the optical field associated with the |c 3 ) quantum state described in the previous 
sections. The LCOS SLM was developed by Kent State Liquid Crystal Institute and 
Hana Microdisplay under the Defense Advanced Research Projects Agency (DARPA) 
Terahertz Operational Reachback (THOR) program.20 The LCOS SLM operates on the 
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far field pupil plane pupil plane 

amplitude amplitude phase 




Figure 3. Numerically calculated amplitudes and phases resulting from holographic 
generation of complex field. The calculation includes the effects of pre compensating 
phase holography artifacts and of spatially isolating the relevant diffracted order. 
Several values of the reference wave tilt are shown. 

principle of discrete-element optical path modulation over a 768 x 1024 element XGA 
format. The optical path length of each of the elements may be independently varied 
over a range of more than 1 micron corresponding to about two waves of modulation at 
the 532 nm wavelength used in this demonstration. Prior to using the LCOS SLM in 
this demonstration, the nonlinearities in the optical phase response were characterized 
and then utilized in a computer algorithm to linearize the phase response to the CGH 
function. In addition, interferometric measurements of the SLM surface flatness were 
utilized to generate a second aberration correction term in the preconditioned CGH 
phase function (p(p,9). The preconditioned CGH phase function is written as, 

(j) ( P7 e) = m-j(kw(p,e)) (40) 

where W(p, 6) is the wavefront error associated with the SLM surface flatness and the 
parameter a' is again taken to be 1.72. The peak-to-valley magnitude of the SLM 
wavefront error is 2.15 waves at 532 nm optical wavelength. With 8-bit resolution, 
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Figure 4. Calculated inner-product-based probability associated with the numerical 
calculation illustrated in Fig. 3. 



the phase response of the SLM spans a range of about 4-ir saturating at a gray scale 
value of 190. The approximately 3.44 radian phase range used in this demonstration 
is accomplished over about 18 phase levels. Details of these procedures have been 
published previously. [20], [21] 

The experimental setup is shown in Fig. 5. The output of a 532 nm continuous- 
wave frequency-doubled Nd:YAG laser is spatially filtered, collimated to approximately 
50 mm beam diameter and propagated to the LCOS SLM. The aperture stop defines 
the area of the SLM utilized to be 1 cm or approximately 90% of the active area. The 
reflected/diffracted optical field is directed to a one- meter focal length lens to create a 
far- field plane where another aperture is inserted to spatially select the m — — 1 order 
of the hologram. The transmitted order is then propagated to a 250mm collimating 
lens that subsequently creates a re-imaged pupil plane where a camera measures the 
pupil plane intensity distribution. The optical phase is measured interferometrically by 
inserting the reference mirror as shown to create a planar reference wavefront and form 
an interferogram on the camera. 

The CGH function of Eq. 20 is calculated using 100 waves of reference wave tilt and 
using the phase function <p(p, 9) in Eq. 40 and the amplitude function a(p, 9) satisfying 
Eq. 29 with A(p,9) and 9) given by the |c 3 ) optical field as given in Eqs. 31. 
The CGH function is scaled to produce a phase range of 3.44 radians on the SLM 
and sampled to define a 768 x 1024 XGA signal that is sent to the SLM driver. The 
experimentally generated optical field may be compared to the theoretical |c 3 ) state as 
follows. The theoretical irradiance function is calculated form the square of the |c 3 ) 
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Figure 5. Schematic of experimental setup showing locations of spatial light 
modulator (SLM), aperture stop (AS), beam splitter (BS), reference mirror (RM), 
objective lens (LI), spatial filter (SF), field lens (L2), and CCD camera. 
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Figure 6. Demonstration of state generation utilizing CGH and LCOS SLM 
technology comparing calculated irradiance and interfcrograms to those generated 
experimentally. 



state amplitude in Eqs. 31. Interferograms associated with the theoretical |c 3 ) optical 
field are calculated by summing the complex field of Eqs. 31 with a planar reference 
wavefront and finding the squared modulus. This is done for 10 waves of interferometer 
reference tilt in both the vertical and horizontal dimensions as shown in the first row 
of Fig. 6. The second row of Fig. 6 shows the experimentally measured irradiance and 
interferograms corresponding to the theoretical case in the first row. The experimental 
results are in qualitative agreement with the analysis presented in this paper. 
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6. Conclusions 

Computer-generated holography with spatial light modulator technology is evaluated 
for its utility in generating QKD states based on photon orbital angular momentum. 
Utilizing azimuthal phase functions to define the orthonormal pure states of a 3- 
dimensional basis, the mutually unbiased bases of a 4-bases system are calculated and 
shown to be azimuthal amplitude and phase functions of the optical field. An analysis 
of CGH with 2-dimensional phase modulators shows that these optical fields may be 
holographically generated and that artifacts associated with thin phase holography may 
be compensated in the computation of the hologram. Furthermore, a numerical analysis 
of practical effects associated with spatially isolating the holographically generated field 
defines a regime in which the azimuthal field may be constructed with high fidelity. 
A demonstration with a LCOS SLM shows practical implementation of this approach 
to generating quantum states for QKD with commercially available technology. Some 
effects associated with SLM technology, such as pixilation and discretization of phase 
levels, were not included in the analysis presented here. These details are important in 
the optimization of SLM designs for implementation of OAM-based QKD. 
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